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Abstract. While many classical traffic models treat the spatial extension of streets 
continuously or by discretization into cells of a certain length Ax, we will subdivide 
roads into comparatively long homogeneous road sections of constant capacity with 
an inhomogeneity at the end. The related model is simple and numerically efficient. It 
is inspired by models of dynamic queueing networks and takes into account essential 
features of traffic flows. Instead of treating single vehicles or velocity profiles, it focusses 
on flows at specific cross sections and average travel times of vehicles. 
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In the last decade, physicists have made many significant contributions to traffic 
modeling pQ. The modeling approaches reach from cellular automata [2 3 over microscopic 
car-following models [H| up to macroscopic or fluid-dynamic traffic models jUEHEj. They 
have, for example, been the basis for understanding the various breakdown phenomena 
of traffic flow in terms of non-equilibrium phase transitions [Zj . 

Classical traffic models are mostly based on the treatment of interacting vehicles, 
their statistical distribution, or their density and average velocity as a function of space 
and time. For numerical reasons, space and time are often discretized into small intervals 
Ax and At. Depending on the model, typical values of Ax are between 5 and 500 meters. 
Although drivers are normally interested in travel times, most models focus on the 
velocities (from which the travel times are sometimes derived). An alternative approach 
are queueing models (see [HI El GDI E] arid references therein): For example, Cremer and 
Landenfeld have developed a microscopic model for saturated urban road networks which 
is discrete in time. It distinguishes only vehicles moving with the maximum allowed 
speed V® = 50km/h and standing vehicles with the space requirement 6.7m, assuming 
discharge rates of Q ont = l/2s and parallel forward motion of all queued vehicles when 
a site becomes empty [H]. The model by Eissfeldt, Grafe, and Wagner eliminates vehicle 
motion and translates macroscopic flows into microscopic time headways. It is based on 
a section travel time and additional traffic-state dependent waiting times of cars, taking 
into account the finite velocity of the upstream propagation of emptied road space |10j . 
It can reproduce empirical flows quite well [TO; • In contrast to these microscopic queueing 
models, Kerner proposes a macroscopic approach to city traffic relating to the continuity 
equation [TT]. It assumes a constant density of saturated (queued) traffic and a time- 
dependent, but space-independent density of unsaturated traffic. Therefore, it neglects 
the spatio-temporal propagation of density variations. 

In the following, we will propose a model which exactly integrates the Lighthill- 
Whitham model over finite road sections. Although it formally eliminates the traffic 
dynamics on the road sections by expressing it through the dynamics at a few cross 
sections of the road network, it takes into account continuous (and discontinuous) 
changes in the density of free or congested traffic and their characteristic finite 
propagation velocities. The model was inspired by a dynamic model of supply and 
production networks [12]. When we now specify road traffic as a queueing system, we will 
take into account essential traffic characteristics such as the flow-density relation or the 
properties of extended congestion patterns at bottlenecks. In fact, traffic congestion is 
usually triggered by spatial inhomogeneities of the road network [7j, and queueing effects 
are normally not observed along sections of low capacity, but upstream of the beginning 
of a bottleneck. Therefore, we will subdivide roads into sections % of homogeneous 
capacity and length Lj, which start at place Xi and end with some kind of inhomogeneity 
(i.e. an increase or decrease of capacity) at place Xi + \ = Xi + L^. In other words, the end 
of a road section i is, for example, determined by the location of an on- or off-ramp, a 
change in the number of lanes, or the beginning or end of a gradient. 
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We will derive our model from the continuity equation 

+ gg^lj = source Terms (1) 
ot ox 

describing the conservation of the number of vehicles. Here, p(x, t) denotes the vehicle 
density per lane at place x and time t, and Qi(x, t) the traffic flow per lane. The source 
terms originate from ramp flows QY^it), which enter the road at place Let us 

define the arrival rate at the upstream end of road section i by Qf r {t) = IiQ^Xi + dx, t), 
where dx is a differential space interval and the number of lanes of road section i. 
Analogously, the departure rate from the downstream end of this section is defined by 
Qf cp (t) = IiQi(xi + i — dx, t). The conservation of the number of vehicles implies that the 
departure rate plus the ramp flow determine the arrival rate in the next downstream 
section i + 

Q7m = QT v {t) + Q? mv {t). (2) 

In order to guarantee non-negative flows, we will demand for the ramp flows that the 
consistency condition — Qf° v (t) < QT mp (t) < Qf^^t) is always met. 

Integrating the continuity equation over x with Xi < x < x-i + \ provides a 
conservation equation for the number iVj(t) = J^ +1 dx Iip(x, t) of vehicles in road section 
i. It changes according to 

= QTV) - Qt cp (t) = Q&(t) + QT^if) - Qt p (t) ■ (3) 

In the following, we will try to express the traffic dynamics and the travel times only 
through the flows at the cross sections Xi. For this, we will take into account the following 
simplified features of traffic flow: (i) For free flow, i.e. below some critical vehicle density 
p cr per lane, the relation between the traffic flow Qi per lane and the density p per 
lane can be approximated by an increasing linear relationship, while above it, a falling 
linear relationship is consistent with congested flow-density data (in particularly, if the 
average time gap T is treated as a time-dependent, fluctuating variable) [TH]- This implies 
Qi(x,t) ~ Qi(p(x,t)) with 

^ l{P> I QT ng (p) = (1 " P//W)/T otherwise. W 
Here, denotes the average free velocity, T the average time gap, and Pjam 

the density 

per lane inside of traffic jams. Moreover, we define the free and congested densities by 

P f rm = Qi/V? and pT*{Q i ) = 0--TQi)p im . (5) 

The quantity Q out = (1 — p cr / Pja,m)/T corresponds to the outflow per lane from 
congested traffic ^3]- Depending on the parameter specification, the model describes 
a continuous flow-density relation (for p cr V^ = Q Q ut) or a capacity drop at the 
critical density p CT and high-flow states immediately before (if p CT V® > Q ou t)- (h) 
According to shock wave theory density variations at place x propagate with velocity 
C(t) = [Qi(x + dx,t) — Qi(x — dx,t)]/[p(x + dx,t) — p{x — dx, t)]. Accordingly, the 
propagation velocity is C = V® in free traffic, and C = — c = — f/(Tp max ) in congested 
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traffic. Therefore, it takes the time period T/ ree = Li/V® for a perturbation to travel 
through free traffic, while it takes the time period Tj Cong = Lj/c, when the entire road 
section i is congested, (iii) Now, remember that congestion in section i starts to form 
upstream of a bottleneck, i.e. at place x i+ i. Let Zj(t) denote the length of the congested 
area and x(t) = Xi + \ — kit) = Xi + Li — the location of its upstream front. Then, we 
have free traffic between Xi and Xi + L { — kit), i.e. Qi(x — dx, t) = Qf r {t — (x — Xi)/V®) 
(considering dx — > 0), and congested traffic downstream of x(t), i.e. Qi(x + dx,t) = 
Qt v {t - (x i+1 - x)/c). With dx/dt = -dk/dt = C(t) and Eq. © we find 

dk = Qt ep (t - li(t)/c)/Ii - QTit - [Lj - k(t)]/V?)/h m 

dt ~ pr g (g dep (i _ k{t)/c)/Ii) _ p fee {Qr{t _ ^ _ l^/V?) / h) ' 

(iv) The capacity of a congested road section i is approximated as the outflow Q out = 
(1 — Per/ P)am)/T from congested traffic per lane times the number Ii of lanes, minus 
the maximum bottleneck strength at the end of this section. This may be given by 
an on-ramp flow Ql amp (t) > or analogously by (Ii — Ii+i)Q ut in case of a reduction 
Ii + i — Ii < in the number of lanes, or in general by some potentially time-dependent 
value AQi(t) in case of another bottleneck such as a gradient: 

QT P (t) = hQ out - max[Qr mP W, {h ~ 4n)Qout, AQ^t), 0] . (7) 

Analogously, the maximum capacity Qf ax (t) of the road section i under free flow 
conditions is given by the maximum flow Iip CT V^ minus the reduction by bottleneck 
effects: 

Qf^it) = I lPcT V? - max[Qr P (t), {U ~ AQ 4 (t), 0] . (8) 

Note that an off-ramp flow Qi amp (t) < does not contribute to these formulas, but one 
may reflect the bottleneck effect of weaving flows (i.e. frequent lane changes at exits) by 
some AQi(t) > 0. (v) We can distinguish three states Si(t) of road section i: Si(t) = 
corresponds to free traffic, which assumes that the length li(t) of the congested area of 
road section i is zero and that the maximum capacity Q^ ax {t — dt) in the last time step 
t — dt, where dt denotes a small time interval, was not reached by the arrival rate at 
time t-dt- 7^ free , i.e. Qf v {t - dt - Tj 1 ' 00 ) < Qf**{t - dt). Si{t) = 1 corresponds to a 
completely congested road section i, which assumes with li(t) = Li that the congested 
area expands over the full section length Li and that the arrival rate Qf r {t — dt) in 
the last time step t — dt was not below the departure rate at time t — dt — Ti° ng , i.e. 
Qf ep (t-dt-Ti° ng ) < Qf x (t-dt). Otherwise we have partially congested traffic in road 
section i and set Si(t) = 2. Altogether this implies: 

f if k(t) = and Qf x {t - dt - T* ree ) < Qf &x (t - dt), 
S l (t) = l 1 if k(t)=Li and Q* cp (t-dt-Tr g )<Q? T (t-dt), (9) 
i 2 otherwise. 

(vi) Finally, we have to specify the departure rate Qf ep {t) as a function of the respective 
traffic situation. Focussing on the cross section at location Xi + \ and considering the 
directions of information flow (i.e. the propagation direction of density variations), we 
can distinguish three different cases: 
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1. If we have free traffic in the upstream section i and free or partially congested 
traffic in the downstream section i + 1, density variations propagate downstream 
and the departure rate Qf ep (t) at time t is given as the arrival rate Qf rr (t — T/ ree ) = 
Qi^Fi(t — T- me ) + Ql^ p (t — T/ ree ), since the vehicles entering section i at time t — T^ TCC 
leave the section after an average travel time T, of T/ ree . 

2. In the case of partially or completely congested traffic upstream and free or partially 
congested traffic downstream, the departure rate Qf ep (t) is given by the capacity 
Q^ p (t) of the congested road section i. 

3. In the case of congested traffic on the entire downstream road section i + 1, the 
departure rate Qf ep (t) is given by the departure rate Qi+i(t — Tj C °° g ) from the 
downstream section at time t — T i cong minus the ramp flow Q|; amp (£) entering at 
location x i+ \. 

Summarizing this, we have 

[ Qf*{t - Tt ee ) if S i+1 (t) ? 1 and $(t) = 0, 

QfP(t) = Q-P(t) if Si+1 (t) ? 1 and S t (t) ? 0, (10) 

{ Qtl\(t - T^) - QT mp (t) if S i+1 (t) = 1. 

A numerical solution of the above defined section-based queueing-theoretical model 
is carried out as follows: First, determine the states of the road sections according to 
Eq. (JHJ). Second, calculate the new arrival and departure rates by means of Eqs. (0) 
and taking into account the boundary conditions for the flows at the open ends of 
the road network. Third, if Si(t) = 0, then set Uit + dt) = l^t) = 0. If Si(t) = 1, 
then set U(t + dt) = k(t) = Lj. Otherwise, if Si(t) = 2, determine the length 
li(t + dt) = kit) + dt dli(t)/dt of the congested area in road section i with formula ©• 
Next, continue with the first step for the new time point t + dt. It is obvious, that this 
numerical solution is significantly more simple and robust than the numerical solution 
of the Lighthill-Whitham model, as shock waves (i.e. the interfaces between free and 
congested traffic) are treated analytically. 

Let us now derive a general relationship for the average travel time Tj(t) of a vehicle 
that enters road section i at time t. The average travel time T,(t) is given by the fact 
that this vehicle will leave section i, when (on average) the iVj(t) vehicles that are in 
section i at time t have passed the downstream end of that section. Since the number of 
vehicles passing the downstream end of section i is given as the time-integral over the 
departure rate Qf cp (t), we have the implicit relationship 

t+Ti{t) t+Ti(t) t 

Ni(t)= J dt>Qt»(t>)= J dt'Qf p (t')- Jdt'Qf^t'). (11) 

t — oo — oo 

Identifying Eq. © with the time-derivative of this equation finally leads to the delay- 
differential equation 

dm = grit) _ = Qt p i(t)+Qn p (t) _ (12) 

dt Qf*(t + T t {t)) ~ Qf p (t + T,(t)) ' 1 ] 
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According to this, the travel time Ti(t) increases with time, when the arrival rate Qf rr 
at the time t of entry exceeds the departure rate Qf ep at the leaving time t + Ti(t), 
while it decreases when it is lower. It is remarkable that this formula does not explicitly 
depend on the velocities on the road section, but only on the arrival and departure 
rates. The calculation of the travel time based on the velocity v(t) is considerably 
more complicated as it depends on the traffic density p(x(t),t) at its respective location 
x(t)=x l + J t t o dt'v(t'). 

In summary, we have proposed a novel queueing theoretical model, which facilitates 
to simulate the departure rates Qf ep (t) of road sections % and the average travel times 
Ti(t) of vehicles in an efficient way. Assuming a fundamental diagram with linear free and 
congested branches (i.e. constant propagation velocities and — c in free and congested 
traffic, respectively) allowed to eliminate the traffic dynamics within homogeneous road 
sections apart from the location of shock fronts, i.e. the moving interfaces between free 
and congested traffic. By determining the travel times from the arrival and departure 
rates of vehicles, the model considerably differs from classical traffic models, in which 
they are determined from the spatio-temporal vehicle speeds. However, note that there is 
some relationship with cell-based approaches [TH] . if we subdivide our large homogeneous 
road sections into smaller subsections (cells). 

Apart from numerical efficiency, the proposed section-based model is promising for 
analytical investigations of traffic in road networks and of dynamic assignment problems. 
Moreover, it can describe the hysteretic breakdown of traffic flow and reproduce typical 
congestion patterns. Nevertheless, it is not as accurate as other macroscopic traffic 
models jSj, which do not assume constant propagation velocities and can describe 
emergent stop-and-go waves. Despite of this, the queueing-theoretical traffic model is 
expected to provide reasonable estimates of the average travel times. 

The above model can be extended in several ways: (i) While the model can describe 
boundary-induced stop-and-go waves (characterized by the temporary disappearance of 
a bottleneck UQontif) ~QT P (t)), emergent stop-and-go waves require the formulation of 
a generalized model containing a dynamic velocity equation, which will reduce analytical 
tractability. (ii) Accidents may be treated by a splitting of sections at the location and 
for the duration of an accident, (iii) Stochastic effects in traffic flows could be easily 
incorporated, (iv) Forecasts of travel times would be possible, if the above model would 
be combined with a model for the prediction of origin-destination flows. In such kinds 
of simulations, it makes sense to extend the model to multi-destination flows [15J. (v) 
Finally, the above model may be generalized to the treatment of urban road networks 
with traffic lights, if we specify the capacities QT v {t) m a time- dependent way which 
reflects the effects of amber and red lights. Note that the actual section capacity does 
not change immediately with the traffic light, but with some delay due to reaction and 
acceleration or deceleration times. 
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